Angle dependent conductance in graphene 
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In this paper, we study a theoretical method to calculate the conductance across a square barrier 
potential in monolayer graphene. We have obtained an analytical expression for the transmission 
coefficient across a potential barrier for monolayer graphene. Using the transmission coefficient 
obtained we have an analytical expression for the conductance. This expression will be used to 
calculate the conductance in the case in which there is a potential barrier, which in our case will 
modelise the behaviour of a top gate voltage of a field effect transistor. Once this analysis has been 
performed we study the scenario in which carriers scatter with the potential barrier with different 
incidence angles and we have found that for any incident angle an effective gap is induced. 



I. INTRODUCTION 

Monolayer graphene is a one atom thick carbon layer 
in a hexagonal honeycomb lattice. From the lattice prop- 
erties one can deduce the band structure and, hence, the 
energy spectrum. 

In graphene, charge carriers are described by the Dirac 
equation instead of Schrodinger's as it is the case of tradi- 
tional semiconductors. This is a direct consequence from 
the fact that graphene has two equivalent triangular sub- 
lattices A and B[l . Therefore, the dispersion relation- 
ship is linear and is expressed hy E = ^.hvpk [2] where 
the positive (negative) sign describes electrons (holes) as 
carriers. Thus, it is deduced that carriers in graphene 
have the same dispersion relationship than massless par- 
ticles with velocity vp ^ c/300. 

According to previous studies, graphene might be 
the ideal material for electronic devices due to 
several properties which are present in graphene such 
as its high mobility (up to 15000cm^/ys[7 ), the large 
scattering length [8 and also because graphene can stand 
a current density which is six order of magnitude greater 
than in copper [9 . The conductance of PN junctions in 
graphene has been studied previously [19 and different 
configurations of junctions NNN, NPN, PPP and PNP 
have been also studied [20 . Experimental measurements 
have been carried out in order to study the conducting 
properties of PNP structures [21 obtained by the de- 
position of a top gate separated from the graphene by 
an air gap. Quantum oscillations of the conductance in 
graphene have been studied both theoretically [22] and 
experimentally [23l [24] . 

Because of its electrical properties, graphene is an in- 
teresting potential material to develop nanodevices us- 
able in technological applications, such as field effect 
transistors (FET). In this work we have modelised a FET 
based on graphene and analyzed its electrical transport 
capabilities based on an exact analytical solution to the 
Dirac equation. Finally, we study the different behaviour 
of the conductance shown by a FET considering that the 



current does not flow perpendicularly to the top gate. 

II. MODEL 

We start our model considering that the carriers in 
graphene can be considered as massless ultrarelativistic 
Dirac particles, since the tight binding hamiltonian for 
graphene leads to a linear dispersion relationship. There- 
fore, graphene carriers are described by the Dirac eigen- 
value equations: 

VFcr-p^(r) =E^(y) (1) 

in where a = {(Jx^cFy) are the Pauli matrices, and ^ is a 
two component spinor. 

Considering an external potencial V(x,y), the hamilto- 
nian is: 

H = vf{(JxPx + cFyPy) + V{x, y) (2) 

which yields 

" ^ ^^^^ ~ ^ " 

In the particular case of V a constant, eventually zero, 
the solution will be a plane wave spinor in the form of 

= 71'*' ( 4 ) 

where 

6k = arctan (5) 

Where the positive (negative) sign will be due to elec- 
tronlike (holelike) regime. 
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FIG. 1. Scheme of the device studied in this paper 



Fig. [T] shows an scheme of the device under consider- 
ation: a rectangular single layer of graphene (blue), with 
two gold contacts which will be the source and drain. 
The device also has a back gate and a top gate which 
is on dielectric layer (Si02, PMMA resist or RX resist, 
for example). The back gate controls the charge car- 
rier density of the sample and the top gate modulates 
the current which flows through the device from source 
to drain passing through a square potential barrier. We 
will suppose that the sample is big enough not to be 
considered as a strip, therefore there will not be edge 
effects. Furthermore, at temperatures close to K, the 
charge carrier density n is proportional to Fermi energy 
squared [27 (n ex E^) and will show a linear dependence 
with the back gate voltage (n ex Vbg) [20 . 

We apply the expressions obtained earlier (Eq. [4| in 
order to calculate the transmission coefficient through a 
square barrier potential Vq and width D^. Each charge 
carrier will have an energy E and the angle of incidence 
will be -7r/2 < (f) < 7r/2. 

The wavefunction will be a linear combination of wave- 
functions in regions I and II and III. 
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FIG. 2. Top panel: Carrier energy and barrier potential de- 
pendence with the position in the sample. Bottom panel: 
scheme of the graphene layer and angles of incidence at the 
barrier potential 



For further discussion, we consider that in zone I the 
amplitude of the incident wavefunction is 1, whereas the 
amplitude of the reflected wavefunction is r. In zone II 
the amplitude of the incident wavefunction is a and the 
amplitude of the reflected wavefunction is b. Finally, in 
zone III only transmitted wavefunction with and ampli- 
tude of value t is to be observed. We define s = sgn{E) 
and = sgn{E — Vq) which are related to the nature of 
the charge carriers (hole or electron regime). 

For convenience, we define the following quantities: 
(/) = arctan , kx = kp cos (/), ky = kp sin (/), 



^, e = arctan [^), = ^ {^)" - 



hv 



In order to obtain the parameters of the wavefunctions 
r, a, 6, t we will impose boundary conditions such as the 
wavefunction must be continuous at the borders of the 
barrier potential. Straightforward algebra, whose details 
are given in reference [33 , yields the expression of the 
transmission coefficient 
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(7) 



As we can see, the transmission coefficient depends upon 
the energy of the charge carrier the angle of incidence 
(j) and the parameters of the barrier potential Vq and D 
whereas it does not depend upon the dimensions of the 
graphene sheet W and L. Also, the transmission coeffi- 
cient is symmetric with respect to the angle of incidence 
(j). According to this expression, in the case of normal 
incidence, the transmission coefficient is the unity with 
independance upon any other parameters. That is the 
case of the Klein paradox for Dirac particles [30] , [31] , [32] . 

The analytical expression obtained for the transmis- 
sion coefficient [33 can be used in order to calculate the 
conductance across a potential barrier. The momentum 
of the charge carriers along the graphene is k = kxi-\-kyj. 
The contribution of each carrier to current density along 
the X direction is given by: 

= T {E,Vo,D,4>)vF^ = 



LW 



-s—TiE,Vo,D,<j>)vFCOs<j> (8) 



The total current is caused by the charges with all 
possible values of and ky. The number of carriers 
with momentum between (kx^ky) and {kx + dkx^ky + 
dky) depends on the number of the available states, and 
its occupation degree, given the Fermi-Dirac distribution 
/o(^,Ai) = (l + exp(^-^)/'=-^)-i. 



The area element in the momentum space is 



LWdk^dky 
47r2 



and the number of electrons in this area is (2 due to the 
spin and 2 due to the duplicity of the valleys of Dirac) 



dueiec = 2x2xLW^^^^/o(£;,/i) = ^^^fo{E,ii)dkxdky 

(9) 

Thus, the current density can be calculated by means of 



If we consider that the current is by both, electrons 
and holes, = — oo and E2 = oo 
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/oo 
T{E,Vo,D,4>)fo{E,fi)EdE (13) 
-CO 

We consider now that the contacts are actually noth- 
ing but ideal reservoirs which stablish equilibrium of the 
distribution of electrons to a certain chemical potential 
/JjL (source) and (drain) yielding two electron fluxes. 
When these two currents are considered, being symmet- 
ric the transmission coefficient as it is, we obtain: 
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(14) 



The later expression can be reduced in some cases : 

• If the bias that exists is large at low temperature, 
the contribution due to the reservoir with lesser 
chemical potential (right) is negligible. It is pos- 
sible to obtain the total density of current if we 
consider only the current that provide the reservoir 
on the left. 

• At low temperature {ksT < Ep — /i), when the 
electrons are hightly degenerated, the Fermi func- 
tion can be approximated by the step function, and 
therefore 
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djx = fxdneiec = -s^^T{E,Vo,D,(l))fo{E,/i) cos(j)dkxdky 

(10) 

We proceed now to use polar coordinates, and use the 
relation E = shvpk 
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Thus, the total current density is 
4e 
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• If the bias is very small, we can approximate 
the Fermi distribution by its Taylor approximation 
around the mean value of the chemical potentials 
/i = ^^^^^ and considering that the difference be- 
tween the chemical potentials is established in fact 
by means of the bias /iL — Mi? = e {Vr — Vl)- 

dfo 

= -e [Vr - Vl){IQ) 

and therefore 
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The minus sign describes the correct direction of 
the electrical current carried for both electrons and 
holes. Therefore, since the conductance is given 
by G = ^ and the intensity is / = jxdy the 
conductance is given by G = -^W. We define the 
effective conductance as Gg// = G/W^ and 
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• At ultra low temperatures (at zero temperature /i ■ 
Ef) we can approximate 



dE 
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so that we obtain 
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In units of e'^/h: 



T{Ef,Vo,D,(P) cos <pd(l) (19) 
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(20) 



We point out that if we express the lengths in nm and 
the energies in meV^ it results hvp = 658.2 meV nm and 
the effective conductance is given in units of /h. 

The conductance will depend on the potential applied 
by the top gate (through the transmission coefficient), 
the potential applied by the back gate (through Ep) and 
the width of the graphene sheet, but not upon its length. 

Furthermore, we observe in (Eqj20| that the maximun 
value of the effective conductance is at T(£^i?, Vb, 1^, ^) = 
1, and therefore 
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By introducing the expresion in (Eqj7|) into (Eqj20| we 
obtain the effective conductance for one square barrier 
potential in graphene: 
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As it is seen, the effective conductance of a graphene- 
based FET depends upon the Fermi energy and the prop- 
erties of the potential barrier created by the top gate 
(height Vq and width D) . 

Figure |3] shows conductance versus Fermi energy. It is 
shown that the curve of the effective conductance has a 
local maximum when Ef = 0.5Vo. At this point, the sign 
of the quantity {Ef — Vq)^ 



El sin^ (f) is changed from 



positive to negative. Therefore, in (Eq 22) the sine of 
this quantity will become a hyperbolic sine, and the con- 
ductance will become smaller. Thus, the resistance will 
increase. This situation will persist until a local mini- 
mum is reached, at Ef — Vq. When Ef is greater than 
Vb, since the transmission coefficient tends to unity, the 
effective conductance grows proportionally to the Fermi 
energy, similarly to (Eq|2T]). We observe oscillations of 
the conductance at Fermi energies which verify Ef < Vq. 
For Ef < 0.5Vb there are several oscillations, whereas for 
values of the Fermi energy such that 0.5Vb < Ef < Vq 
there is only one oscillation. If Ef » Vq the conduc- 
tance does not depend on the value of the barrier width 



and varies linearly with the Fermi energy. 



As we can see in figure [4j when Vb increases, the effec- 
tive conductance decreases linearly and takes the same 
value indepently of the value of the width, for Vq veri- 
fying Vq < Ef- When Vb takes the value of the Fermi 
energy, there is a local minimum, whose value depends 
on the width of the barrier D. For Vb > Ef, the conduc- 
tance increases and exhibits oscillations which are bigger 
as the width of the barrier decreases. Furthermore, it 
is appreciated that, as the value of D decreases the am- 
plitude of the oscillations increases. These oscillations 
appear only for values of Vq which verify Vq > Ef- 



The dependance of the effective conductance upon the 
width D of the barrier is shown in figure [5] the con- 
ductance decreases until it tends to be stabilised in the 
vicinities of a constant value. It is also observed that the 
oscillations tend to be lesser as the Fermi energy grows. 
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FIG. 3. Effective conductance versus Fermi energy at height 
of the barrier Vb = 200 meV and width of the barrier D 
50(red), lOO(green) and 150(blue) nm. Conductance has a 
local maximum at Ep — 0.5Vb, a local minimum at Ep — Vo 
and increases linearly with Ep > Vq. 
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FIG. 4. Conductance versus height of the barrier at Fermi 
energies Ep = 100 meV and D 50(red), lOO(green) and 
150 (blue) nm. The three curves take the same values for 
Ep > Vb, and present a minimum for Ep = Vb. The ampli- 
tude of the oscillations depends on the value of the width of 
the barrier D. 

III. WEIGHTED EFFECTIVE CONDUCTANCE 

Up to now we have considered that the carriers scatter 
through the barrier with the same probability for each 
angle of incidence. It should be considered another sce- 
narios in which, due to the inhomogeneities of the mate- 
rial or the design of the device, the carriers flow around 
certain direction. Thus, the angle of incidence will be 
different from 0° and will obey certain probability distri- 
bution function P(0) centered around an angle 0o- This 
probability distribution function must verify the condi- 
tion 

/ P^dcl) = TT (23) 
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FIG. 5. Conductance versus width of the barrier when assum- 
ing a height of the barrier Vb = 100 meV and Fermi energies 
Ep = 25(red), 50(green), 75(blue), lOO(yellow) and 150(pur- 
ple) meV in order to perform the calculations. 



Next, we will consider different probabilities of distri- 
bution and we will study its effect on the effective con- 
ductance. 



A. Gaussian distribution 

We consider now the case of a Gaussian function for 
the probability distribution of the angle of incidence. 
This assumption represents that transport is of ballistic 
nature and that there are a few inhomogeneities which 
cause scattering processes within the sample and devi- 
ates carriers from their ballistic trajectory; making them 
to impact the barrier potential with a non zero angle of 
incidence. Since not all the carriers will suffer the same 
scattering process, we can assume that the majority of 
carriers will impact the barrier potential with a similar 
angle and that there will be more carriers impacting 
the barrier potential with an angle of incidence 0o than 
carriers impacting the barrier with an angle of incidence 
that differs substantially from 0o- Thus, we assume 
that the angles of incidence of the carriers will obey 
a probability distribution function that has the form 

P(0)=cexp(-(^). 

The parameter a of the Gaussian functions makes the 
distinction between an homogeneus sample (a narrow 
Gaussian bell) or a sample with plenty of inhomo- 
geneities (a wide distribution). A narrow Gaussian 
function represents the case in which almost all the 
carriers scatter with the barrier at the same angle 
of incidence. This will happen when the sample is 
homogeneus and the carriers do not scatter with the 
sample inhomogeneities, and therefore the transport will 
be ballistic and almost all the carriers will move through 
the sample with a similar direction. Obviously, a wide 
Gaussian function represents the opposite case. 
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We start the analysis by considering the scenario in 
which the transmission coefficient equals the unity: <po 
=0°. 

Obviously, the parameters a and c are not indepen- 
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FIG. 6. Effective conductance versus the Fermi energy when 
the distribution function of the angle of incidence is a Gaus- 
sian function centered around 00 = 0° assuming different 
widths of the Gaussian: 0.64 (red), 0.21 (green) and 0.08 
(blue). The height of the barrier is Vb = 200 meV and its 
width is D — 100 nm. 



dent, since they must hold the normalization condition 
acy^Erf (^^^^^ = TT. As we can see, as c becomes 
smaller a becomes greater. 

The conductance is given by: 



(24) 
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FIG. 7. Conductance versus the height of the barrier Vb 
asuming a Gaussian probability distribution function cen- 
tered around 0o = 0°. The width of the Gaussian distri- 
bution is a = 0.64 (red), a = 0.21 (green) and a = 0.08 
(blue). The Fermi energy of the system is assumed to be for 
Ef — 100 meV and the width of the barrier is D = 100 nm. 
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As in the previous case, the effective conductance 
depends on Ep^ Vq and D but it shows a different 
behaviour, as shown in figures |6] and jT] 

We study the behaviour of the conductance when the 
Fermi energy is changed while the other variables remain 
at fixed values (see figure [6]). The conductance oscillates 
whilst it increases until a local maximum is reached. 
The value and position of this maximum depends on 
the parameters of the Gaussian function used. Once 
the maximum is reached, the conductance decreases 
until it reaches a local minimum when Ep ~ Vo- Once 
the minimum is surpassed, the conductance increases 
monotonically. 

At it is shown in figure [7) when using Vb as the driving 
parameter and the other parameters remain fixed at con- 
stant values, we can see that the effective conductance 
decreases until it reaches a minimum around Ep ^ Vq. 
It is also observed that the narrower Gaussian function 
(which is equivalent to a less disordered distribution) 
the variation of the conductance is more pronounced. 



Once Ep \s greater than Vb we can observe oscillations 
of the conductance whose frequency tends to the value 
uj = 2D/hvF. 

We consider now the case in which the Gaussian dis- 
tribution is not centered around (j)o = 0°. This might 
be done by tilting the orientation of the top gate in a 
graphene based device as, for example, shown in [29] . 
In this case, the normalization condition is given by the 

expression ac^ (Ev/ (^ V+ Erf (^)) = n. 

As we can see in figures [9| and \lO\ the results are similar 
than the obtained in the previous case, but the minima 
observed previously are wider. In fact as a is closer to 
zero these minima become wider. 



B. Delta distribution 

Next, we consider the extreme scenario in which the 
Gaussian function is so narrow that it can be represented 
by a Delta function S{(j) — This is the case in which 
the material is ideal and all the carriers scatter through 
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FIG. 8. Conductance versus the width of the barrier using 
a Gaussian distribution of probabihty centered around 00 = 
0°. The width of the Gaussian is assumed to be a = 0.64 
(red), a = 0.21 (green) and a = 0.08 (blue). The Fermi energy 
is Ef = 50 meV (hne) and Ef = 150 meV(dashed Une). Also, 
the height of the barrier is Vb = 100 meV. 
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FIG. 9. Conductance versus the Fermi energy considering 
that the distribution is a Gaussian function centered around 
00 = 22.5° using different widths for the Gaussian function: 
a = 0.65 (red), a = 0.21 (green), a = 0.08 (blue). The height 
of the barrier is Vb = 200 meV and its width D — 100 nm. 



the barrier with the same angle of incidence 0o- 
In this case, the conductance is given by: 

Grf,0o(^F,^o,l^) = ^|^T(Ef, ^0, A 0o) COS 00 (25) 

It is observed that the effective conductance depends 
now on the angle of incidence, along with the rest of the 
parameters previously considered (Fermi energy, Ep^ 
height of the barrier, Vb, and width of the barrier, D). 

In particular, if the current flows perpendicularly to 
the barrier, the transmission coefficient equals 1 with 
no dependance on the Fermi energy and the param- 
eters of the barrier, as a consequence of the Klein 



FIG. 10. Effective conductance versus the height of the bar- 
rier Vb considering that the distribution is a Gaussian func- 
tion centered around 0o = 22.5° using different widths for the 
Gaussian function: a = 0.65 (red), a = 0.21 (green), a = 0.08 
(blue). The Fermi energy is Ef — 100 meV and its width 
D — 100 nm. We note that as the Gaussian becomes nar- 
rower, i.e, as the number of inhomogeneities is decreased, the 
minimum of conductance is wider. 
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FIG. 11. Effective conductance versus the width of the barrier 
D assuming that the distribution is a Gaussian function cen- 
tered around 0o = 22.5° for a = 0.65 (red), a = 0.21 (green), 
a = 0.08 (blue) for Ef = 50(line), 150 meV(dashed hne) and 
= 100 meV. 



paradox [30] , j3T] , [32] . Therefore, the conductance is pro- 
portional to the Fermi energy: 



\Ef\ 



For non normal incidence, the conductance presents an 
effective gap between the values of the Fermi 

energy, developing an effective minimum at Ep = 
Vo/ cos^ 00 (as shown in figure 12), whose value is 
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In this scenario in which we assume a delta-like proba- 
bility distribution function, we can see that for an angle 
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FIG. 12. Effective conductance versus Fermi energy assum- 
ing a delta like distribution when the barrier is assumed to 
be such that its height is Vb = 200 meV and its width 
D = 100 nm. We assume different angles of incidence 00 = 
0°(red), 2° (green) and 5° (blue). These curves are the limit of 
a narrow gaussian. 




0.05 - 



0.00 F - 

50 100 150 200 250 300 




0.1 - 



O.Q l ^ 

50 100 150 200 

D (mil) 

FIG. 14. Effective conductance versus the width of the barrier 
if a delta like distribution probability is assumed at Fermi 
energies Ep = 50(line), 150 meV(dashed line), and the height 
of the barrier potential is Vb = 100 meV. We have assumed 
different angles of incidence: 00 = 0°(red), 2° (green) and 
5° (blue). 
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FIG. 13. Conductance versus the height of the barrier assum- 
ing a delta like distribution supposing that the Fermi energy 
is Ef — 100 meV and that the width of the barrier \s D — 
100 nm. This calculation has been performed assuming differ- 
ent angles of incidence 00 = 0°(red), 2° (green) and 5° (blue). 



FIG. 15. Effective conductance versus the Fermi energy when 
a delta-like probability distribution function is assumed at 
angle of incidence 0o = 22.5° and when the potential barrier 
is such that Vq = 200 meV and D = 100 nm. 



of incidence of 0o = 22.5° (or 7r/8) a wide minimum of 
the conductance is developed (figures [T2| and 13). This 
minimum is formed both when using the Fermi energy 
(which can be associated with the back gate voltage of 
the transistor shown in figure [T]) and the height of the 
barrier (the top gate voltage in our model schematized) 
as a driving parameter. In this case, the dependance of 
the conductance on the width of the barrier shows no 
remarkable features but a small oscillation whose ampli- 
tude depends on the Fermi energy (see figures 14 and 

Trl). 



As the angle of incidence is increased, the features de- 
scribed tend to be greater. In particular, the effective gap 
which appears is wider and more pronounced, as seen in 



incidence 0o = 22.5°. The observed minimum of the 
conductance shows that in certain conditions the Klein 
tunneling might be avoided and that the fabrication of 
an effective FET based on graphene with a noticeable 
on-off ratio is feasible. 



IV. CONCLUSION 

We have obtained an analytical expression of the 
transmission coefficient through a square barrier poten- 
tial based on the previous work done [33] in which the 
continuity conditions were applied in order to obtain the 
Dirac functions of carriers in graphene in the presence 
of a barrier potential. 



figures 15 and 16 where we have assumed an angle of 
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FIG. 16. Conductance versus the height of the barrier Vb 
when a delta hke probabihty distribution function is assumed 
at angle of incidence (j)o = 22.5° and the other parameters are 
such that Ef 100 meV and D 100 nm. 
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FIG. 17. Conductance versus the width of the barrier assum- 
ing a delta like function centered around 00 = 22.5° for Ep 
= 50(line), 150 meV(dashed line) and Vo = 100 meV. 



Within this approach, we have stablished a theoretical 
model to study transport through a square barrier po- 
tential in graphene. Analytical expressions for the trans- 
mission coefficient in several different scenarios have been 
obtained. These different scenarios are a modelization of 
the different scattering conditions of the sample due to 
the different degree of inhomogeneity in the sample. 



In this model we have also studied a graphene based 
device and its effective conductance between its two 
terminals. In particular, we have found an scenario in 
which the Klein paradox is neglected and the creation of 
an effective gap is induced. Therefore we think that it 
might be suitable to control the transport in a graphene 
based field effect transistor. 



We believe that a top gate which is non-perpendicular 
to the source-drain direction might create a gap in the 
case in which the graphene is clean and shows almost 
no inhomogeneities. This phenomenon occurs for any 
angle of incidence The position and width of the 
energy gap, depend on both the barrier height and 
itself. This feature might be an important step in the 
development of transistors based on graphene, as the 
non existance of a gap is one of the main difficulties 
found to effectively modulate the current via the bias. 



The authors thank Jose Maria Cervero for useful dis- 
cussion. The authors greatly acknowledge the finan- 
cial support of this research from the DGICYT under 
project FIS2009-07880 and JCYL SA049A10-2 and MEG 
through grant AP2009-2619 



[1] A.H. Castro Neto, F. Guinea, N.M.R. Peres, K.S. 
Novoselov, A.K. Geim, arXiv:0709.1163v2 [cond-mat] 
(2008) 

[2] M.I. Katsnelson, Mater. Today 10, 20 (2007) 
[3] J.M. Cervero and E. Diez, Int. J. Theor. Phys. 50, 2134 
(2011) 

[4] M.C. Lemme, T.J. Echtermeyer, M.Baus, H. Kurz, IEEE 

Electron Device Lett. 28, 282 (2007) 
[5] Y.-M. Lin, C. Dimitrakopoulos, K.A. Jenkins, D.B. 

Farmer, H.Y. Chiu, A. Grill, Ph. Avouris, Science 327, 

662 (2010) 

[6] F. Xia, D.B. Farmer, Y.-M. Lin, Ph. Avouris, Nano. Lett. 
10, 715 (2010) 

[7] A.K. Geim, K.S. Novoselov, Nature Materials 6, 183 
(2007) 

[8] K.S. Novoselov, A.K. Geim, S.V. Morozov, D. Jiang, Y. 
Zhang, S.V. Dubonos, I.V. Grigorieva, A. A. Firsov, Sci- 
ence 306, 666 (2004) 



[9] A.K. Geim, Science 324, 1530 (2009) 
[10] N.M.R. Peres, J. Phys.: Condens. Matter. 21, 323201 
(2009) 

[11] J. Milton Pereira Jr., P. Vasilopoulos, F.M. Peeters, Mi- 
croelectronics J. 39, 534 (2008) 
[12] H.C. Nguyen, V.L. Nguyen, J. Phys.: Condens. Matter. 

21, 045305 (2009) 
[13] M.I. Katsnelson, Eur. Phys. J. B 51, 157 (2006) 
[14] K. Ziegler, Phys. Rev. B 75, 233407 (2007) 
[15] K.S. Novoselov, A.K. Geim, S.V. Morozov, D. Jiang, M.I. 
Katsnelson, I.V. Grigorieva, S.V. Dubonos, A. A. Firsov, 
Nature 438, 197 (2005) 
[16] E.H. Hwang, S. Adam, S. Das Sarma, Phys. Rev. Lett. 

98,186806 (2007) 
[17] K. Nomura, A.H.MacDonald, Phys. Rev. Lett. 98, 

076602 (2007) 
[18] M. Trushin, J. Schlieman, EPL 83, 17001 (2008) 
[19] T. Low, S. Hong, J. Appenzeller, S. Datta, M.S. Lund- 



10 



Strom, IEEE Transactions on Electron Devices 56, 
1292(2009) 

[20] B. Huard, J.A. Sulpizio, N. Stander, K. Todd, B. Yang, 
D. Goldhaber-Gordon, Phys. Rev. Lett. 98, 236803 
(2007) 

[21] R.V. Gorbachev, A.S.Mayorov, A.K. Savchenko,D.W. 

Horsell, F. Guinea, Nanoletters 8, 1995 (2008) 
[22] V.A. Yampolskii, S. Savelev, F. Nori, New J. of Phys. 

10, 053024 (2008) 
[23] A.F. Young, P.Kim, Nature Phys. 5, 222 (2009) 
[24] M. Begharbekov, O. Sul, N. Ai, E.- H. Yang, S. Strauf, 

Appl. Phys. Lett. 97, 122106 (2010) 
[25] L Jung, D.A. Dikin, R.D. Finer, R.S. Ruoff, Nano Lett. 

20, 40 (2008) 



[26] N.R. Malik, Circuitos electronicos. Andlisis, simulacion y 

diseno, Prentice Hall (1996) 
[27] H.-S. Philip Wong, D. Akinwande, Carbon nanotuhe 

and graphene device physics, Cambridge University Press 

(2011) 

[28] S. Dditta., Electronic transport in mesoscopic systems, 

Cambridge University Press (1995) 
[29] T. Low, J. Appenzeller, Phys. Rev. B 80, 155406(2009) 
[30] A. Calogeracos, N. Dombey, Contemp. Phys. Rev. 40, 

313 (1999) 

[31] M.L Katsnelson, K.S. Novoselov, A.K. Geim, Nature 

Phys. 2, 620 (2006) 
[32] C.W.J. Beenakker, Rev. Mod. Phys. 80, 1337 (2008) 
[33] J. Lejarreta, C. Hdez. Fuentevilla, under preparation 



